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Abstract 

We show that fermion charge fractionalization can take place in a recently proposed chiral gauge 
model for graphene even in the absence of Kekule distortion of the graphene honeycomb lattice. 
We extend the model by adding the coupling of fermions to an external magnetic field and show 



o 

that the fermion charge can be fractionalized by means of only gauge potentials. It is shown that 
the chiral fermion charge can also have fractional value. We also relate the fractionalization of the 
fermion charge to the parity anomaly in an extended Quantum Electrodynamics which involves 
vector and vector-axial gauge fields. 
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I. INTRODUCTION 



It is known that the single-particle spectrum for the charge carriers on a honeycomb lattice 
described by means of a tight-binding Hamiltonian contains two zero-energy Dirac points. 
The tight-binding Hamiltonian, linearized around the Dirac points, reveals that the kinetics 
of the charge carriers is governed by a single free Dirac equation for a two-flavor spinor whose 
components designate massless fermions on the two triangular sublattices of the honeycomb 
lattice. A convenient coupling of the Fermi surface charge carriers to distortions of the lattice, 
which preserves the invariance of the system under time and space reversion symmetries, 
opens a mass gap in the single-particle energy dispersion relation. Recently, Hou, Chamon 
and Mudry |l| have shown that when such lattice distortions exhibits a vortex-like profile, 
there are zero modes excitations in the single-particle energy spectrum. The vorticity of the 
distortion, either n > 1 or n < -1 , determines which one of the sublattices supports \n\ 
zero modes. From the existence of fermion zero modes and from the sublattice symmetry, 
they show that the fermion quantum charge is fractionalized, similarly to what was shown 
to occur in one dimensional systems such as polyacetylene 4|. 

Jackiw and Pi Q have proposed a chiral gauge theory for graphene. They show that 
the system considered in l| is invariant under global chiral gauge transformation and the 
invariance of the system under local chiral gauge transformation (Ua(1)) is preserved if the 
space component of a vector-axial gauge potential coupled in a chiral manner to the fermions 
is added to the system. Such a naturally incorporated gauge potential is the vector field 
which enters the phenomenological Nielsen-Olesen-Landau-Ginzburg-Abrikosov model of a 
complex scalar field minimally coupled to the U(l) gauge field potential and whose classical 
minimum energy solutions exhibit vortices profiles as those binding zero-energy fermions in 
the scenario for fermion charge fractionalization in In this way, Jackiw and Pi have 
specified the dynamics for the complex scalar field, which plays the role of the distortion of 
the graphene lattice. Moreover, they show that sublattice symmetry is preserved and that 
the fermion zero modes, as well as the fermion charge fractionalization are persistent even 
when the self-consistent vortex solution for the gauge field is incorporated to the system. 

Approximately two decades ago, Semenoff presented a physical realization of the parity 
anomaly in 2+1 space-time dimensions also in a honeycomb lattice J(J. This time the lattice 
contains two species of atoms and could describe for example a monolayer of boron nitride 
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compound. The difference in energy of electrons on the atoms leads to a mass term for 
the fermions in the effective Dirac Hamiltonian in the continuum limit. The addition of 
an external gauge field coupled minimally to the fermion current leads to an abnormal 
induced current for each one of the fermion species, since the fermion mass breaks the parity 
symmetry. Such abnormal current is proportional to the sign of the mass of the fermion, 
and it is persistent when the mass is taken to zero. Since the sign of the mass of one 
species of fermion is contrary to the other one, the total fermion current is null, but an odd 
combination of the different currents leads to a nonzero abnormal current. 

For each one of the fermion species in the system considered in 6|, the induced charge in 
the presence of an external magnetic field is due to the localized fermion zero modes. The 
induced charge may be fractional and proportional to the magnetic flux on the plane [Z| . In 
the case of a vortex magnetic field one has a finite magnetic flux and the number of fermion 
zero modes is proportional to the magnetic flux |8|. 

n 

We show that a vortex magnetic field, as in the system considered in [5fl, binds, by itself, 
zero-energy fermion states whose spinor structure is equal to that found in the fermion 
zero modes in fl and that is, the zero-energy states are supported on only one of the 
triangular sublattices. The breaking of the sublattice symmetry by the fermion zero modes 
implies into an induced fractional fermion charge which is proportional to the magnetic flux. 
In fact, one can see from the spinor structure of zero modes that they are eigenstates of a 
Dirac matrix which anticommutes with the Dirac Hamiltonian and that, the introduction 
of a staggered chemical potential term proportional to that Dirac matrix reveals the origin 
of the sign of the fractional charge. Whether the chemical potential is positive or negative 
implies into the contribution of positive or negative energy states. In this way the realization 
of induced fractional fermion charge in graphene is provided in a manner very similar to that 
one proposed originally in 0], that is, by means of only vector-axial gauge field. 

Furthermore, by imposing the invariance of the Dirac Hamiltonian under the usual U(l) 
local gauge transformation, we extend the system by minimally coupling the fermions to 
a vector gauge potential. We show that a fractional fermion charge may also be induced 
by a vortex-like magnetic field associated to the added vector potential with and without 
the magnetic field associated to the original vector-axial gauge potential. We analyze the 
discrete symmetries of the isolated fermion zero modes and find that they are supported on 
both triangular sublattices. Based on that symmetry, we show that there is room to the 
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fractionalization of the fermion (chiral) charge as well. 

In the next section we consider the Dirac Hamiltonian resulting from the linearization, 
around the Dirac points, of the tight-binding Hamiltonian for the fermions on the graphene 
honeycomb structure. In this case we start with the extended Dirac Hamiltonian proposed 
in J|]] and in the absence of the complex scalar field and show the isolated fermion modes as 
well as the induced fermion charge in the presence of a vortex magnetic field associated to 
the vector-axial gauge potential alone. Next, we extend the system by coupling the fermions 
to the usual vector gauge potential and show the isolated fermion modes. By analyzing 
the discrete symmetries of the problem we calculate the induced fermion charges in the 
presence of the two kinds of magnetic vortices (or solenoids). The third section is dedicated 
to a detailed analysis of the fermion states when the magnetic fields are homogeneous and, 
from the behavior of the lowest Landau level (LLL), we illustrate the origin of the induced 
fermion charges. In the fifth section we obtain the effective Chern-Simons (C-S) action in 
an extended Quantum Electrodynamics (QED)in 2+1 dimensions with vector and vector- 
axial gauge fields and from that we calculate the associated induced fermion currents which 
recovers the fermion charges obtained in the previous sections. The sixth section is left to 
the conclusions. 

II. FERMION ZERO MODES ON VORTEX-LIKE MAGNETIC FIELDS 

The field theory Hamiltonian describing the dynamics of the electrons on the graphene 
honeycomb structure in the presence of a complex scalar field, which play the role of the 
lattice distortions, together with a vector gauge potential coupled in a chiral way to the 
electrons is written in a simplified version as 



H 



I 



d 2 r^\r) 



(1) 



where ty(f) is the four component spinor 




(2) 
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and a, and 75 are Dirac matrices in the representation explicit below. The superscripts 
in the spinor components designate the triangular sublattice, A or B, the electrons are 
supported on, while the subscripts stand for each one of the Dirac points (the single-particle 
energy spectrum obtained from the original tight binding Hamiltonian exhibits two Dirac 
points). The Kekule distortion is represented by the complex field ip(r) = \ip(f)\e lx ^ = 
(p r (r) + iip 1 if) and A^(r) is the vector- axial gauge potential added to make the Hamiltonian 
invariant under local chiral gauge transformation 

matrix structure in ([1]) is revealed by means of the following gamma matrices 

? 





75 = — ia x a 2 a 3 = I I , (3) 

-/ 

where I is the 2x2 identity matrix and ~a = (a±, a 2 ) and 03 are the Pauli matrices. 

Now, we consider the Hamiltonian ([I]) by dropping out the interaction term of the elec- 
trons with the Kekule distortion and show that if the vector-axial gauge potential has a 
convenient asymptotic behavior, there will be isolated fermion zero modes in the single- 
particle energy spectrum. For the sake of simplicity we consider 

A\ = jtdjMr) ^B 5 = -V 2 A(r), (4) 

where A$(r) is a scalar function of the radial coordinate only. The time-independent Dirac 
equation 

a.(-iV - <?75^Mr) = EV(f), (5) 

can be rewritten as four equations for the spinor components (the coupling constant q > 
is introduced here for future convenience). Two of them, in a simple form, are 

- 2i [d- z - qd- z As] < = mi and - 2i [d z + qd z A 5 ] ^ = Eip b +J (6) 

whilst the equations for i\) h _ and xfit are obtained by taking the complex conjugate of the 
above equations and by making the identifications ip_ = {ip b + )* and %pt = (^+)*- I n the 
above equations we are using the notation z = x + iy and d z = (d x — id y )/2. In such a 
simple notation, one can easily check that the zero-energy non-normalized eigenstates would 
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be given by 



^ = f( z )e qM{r \ ipl = h(z)e- qA5{r \ it = f(z)e qA ^'\ i\) a _ = h(z)e~ qA5{r) , (7) 

where f(z) and h(z) are holomorphic functions. The normalization of the eigenfunctions 
depends on the asymptotic behavior of A^{r) and, from the functional dependence of the 
spinor components, either ^>±(f) or ^±(r) is normalizable, thereby one obtains one of the 
following zero-energy eigenstates 

/ e +ue\ ( o \ 








r l e qA${r) 



or * (r) = N 



-no 



r l e -qA 5 {r) 



(8) 



V o J 



where I is taken to be an integer number and 9 is the angular variable in cylindrical co- 
ordinates. The degeneracy of the zero-energy level depends on the number of values I can 
assume and it is linked to the magnetic flux according to the theorems proved in 8|. In this 
section we are concerned with finite magnetic flux configurations only. 

To show that a fractional fermion charge can be induced with such magnetic field con- 
figurations, we note that the zero modes are eigenfunctions of a 3 and that the Hamiltonian 
operator a.(p — 75A5) anticommutes with the matrix a 3 , then one can show that each 
eigenfunction of negative energy is obtained from an eigenfunction of positive energy by 
\I/_|E|(r) = a 3 ty\E\ In other words, the energy spectrum is symmetric around the zero- 
energy level whose corresponding eigenfunctions are self-conjugate. From this symmetry 
and from the fact that the zero modes contain representatives of only one of the triangular 
sublattices one obtains, by following the procedures in Q, Q and that the induced 
fermion number is given by 



N 



d 2 f -(0|[* f (f),*(f 



T 4vr ' 



(9) 



where $5 = J d 2 rB^ is the magnetic flux and the upper (lower) sign on the right-hand side 
of ([9]) comes out when the zero modes are assigned to fermions (antifermions). 

The induced fermion charge is finite provided the magnetic flux is finite and that can 
be achieved by means of for example a solenoid or vortex magnetic field, such that [3] 



A&(r — > 00) 



-($ 5 /27r) lnr. For $ 5 > one has v&jj (r -)■ 00 



~ r aw (e 



-ue 
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and 
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for $ 5 < (r — > oo) ~ r /+ ^r"(0 e~ lW e +lW 0), where we are assuming that / < [|^| — 1], 
with [v] denoting the largest integer less than v. The charge is fractional if the magnetic flux 
is an integer odd number (n ) of quantum of flux hvp/q, where h is the Planck constant, Vf 
is the Fermi velocity and q is the coupling constant associated to the chiral coupling, that 
is, 

Q = eJ d 2 r i<0|[*t(r),*(f)]|0) = =F« |, (10) 
where e is the unit of electric charge. 



A. Minimal coupling to a vector gauge potential, zero modes and fermion charges 
fractionalization: 

The Hamiltonian ([1]) is also invariant under the usual global gauge transformation \1> — >■ 
e l *ty, with £ uniform and constant. By promoting such a transformation to a local one, a 
vector gauge potential A must be added to the kinetic term in such a way that 

We are now concerned with the appearance of zero-energy bound states of the system 
when such a vector gauge potential is also taken into account, namely we are going to 
consider the Dirac equation 

a.(-iV - rfc^ - eA)^{f) = E^(r), (11) 

particularly when both gauge potentials engender finite magnetic flux. Here we have also 
explicit the coupling constant e > 0. 

An interesting aspect of the Hamiltonian operator in (ITT]) is that it anticommutes with 
a 3 and with a 3 7s. In the previous case the possible zero-energy solutions are eingenstates 
of a 3 , whilst in the case A 5 = and ~X given by 

A { = e ij djA(r) ^B = -V 2 A{r), (12) 

the zero-energy solutions are eingenstates of a 3 75. In this case the spinor components satisfy 
the equations 

- 2i [d z + ed z A] = - ^ [d- z - ed- z A] ^ = 

2i [d z + ed z A] it = Evj)*, 2i [d- z - ed- z A] = Eip b _, (13) 
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whose eigenfunctions for E = are given by 

/ e +ue \ 



Mr) =M 





,+ue 



( \ 

„-iW 



r l e eA(r) 



or V (r) = N 



\ o J 







\ e~ lW ) 



r l e -eA(r) 



(14) 



We notice that fermion zero modes are supported on both sublattices and are eigenstates 
of a 3 75. The energy spectrum is symmetric around E — 0, the negative energy eigenstates 
are obtained from the positive energy ones by means of the norm preserving operation 
^_|E|(r) = a 3 75\l/|E|(r). Thus, the induced fermion number is given as in ([9]) with g$5 
replaced by e$, with $ = J d 2 fB. By considering A(r — > oo) ~ — ($/2-7r)lnr we have 
tfjftr oo) ~ r'-if (e +iW e +i/e 0) for $ > and ^(r oo) ~ r l+ £ (0 e~ iW e" i/e ) 
for $ < 0. An induced fermion charge is obtained similarly to the previous case and if the 
magnetic flux is quantized as an integer odd number of quantum of flux hvp/e the fermion 
charge is fractionalized and given by (II Op . 

The fractionalization of the fermion number (charge) is also possible in very special 
circumstances when both quantized magnetic fluxes acts simultaneously. Moreover, another 
charge - we call it chiral charge - can also be fractionalized. That is what we are going to 
show next. 

The spinor components satisfy the following equations 



- 2% [d z + d z A+] r+ = E^j b + , - 2i [d- z - d- z A + ] ^ b + = ^ 
2% [d z + d z AJ\ ip b _ = Ei) a _, 2i [0 - d^A-} ipl = Eijt, 



(15) 



In the equations above we have set H — Vf — 1, A± = eA ± qA§ and, for sake of simplicity, 
we consider the following asymptotic behaviors A$(r — > oo) ~ — — lnr and A(r — > oo) ~ 
— - lnr, with and n positive integer numbers. 

Normalized zero-energy states are obtained in some particular cases, namely 



/ 



n>n 5 
0,1+1- 



(f) 





1 

y/47Tf™r> 2l - +1 e^-(r>)dr> 





r l+ e il+9 e A+(r) \ 



n > n 5 , l± = 0, 1, [n ± n 5 - 1], 



(16) 



or 






1 



Af(r) \ 



Note that a 3 7 5 ^"> n ^(f) = *of1_(r) and a 3 ^ 



r l- e -il-6 q—A- (r) 



3iT/ n<n 5 



,l+,L 



n < n 5 , l± — 0, 1, [715 ±n — 1]. 



(17) 

^™f n ? (f). On the other 



hand, when the magnetic fluxes are both negative, the zero modes are a 75$ q 



|ra|>|n 5 | . 



l+,L 



1(f) and a*KrX l (r) = < 



n\<\ns 



\n\<\n 5 \ . 



l+,l. 



Both kinds of zero modes given by expressions (Tl6|) and ffTT|) lead to induced fermion 
numbers, which depend on the magnetic fluxes as 



and 



N 



N 



1 n 
d 2 f -(0\[&(r),V(r)]\0) = T~, for | 



(i 2 f-(0|[* t (rl,*(r)]|0} = Ty, for \n\ < \n 5 \ . 



(18) 



(19) 



Again, the upper (lower) sign on the right-hand side of (1T8|) and (Tl9|) depends on the zero- 
energy states are taken as fermions (antifermions). 

We have also found that a fractional (chiral) fermion number can also be induced, namely 



and 



d 2 r -(0|[^(r),75*(r 



n 5 f II II 

T—, for |n| > |n 5 | 



n C 11 11 
=F-, for |n| < 1 7i 5 1 



(20) 



(21) 



if \n\ and jrasj are odd numbers. 

This (chiral) fermion number is proportional to the non-vanishing fractionalf anomalous) 
fermion charge found previously in the Boron-Nitride honeycomb structure 6| where two 
species of electrons which belong to the two different atoms implies into a mass term. Such 
a mass term, in the representation we have been working here, is equivalent to a staggered 
chemical potential. In the next section we analyze in detail the single-particle states for 
electrons in the honeycomb structure and in the presence of homogeneous magnetic fields 
when a staggered chemical potential and a parity-breaking mass term are added to the 
Hamiltonian. 
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III. LANDAU LEVELS AND INDUCED CHARGES ON UNIFORM MAGNETIC 
FIELDS 



In this section we obtain the fermion bound states when the electrons on the honeycomb 
structure are under the action of vector gauge potentials for homogeneous magnetic fields. 
In this analysis we introduce terms in the Hamiltonian density which provide gaps in the 
single-particle energy spectrum. Such terms are a staggered chemical potential fia 3 and 
a parity breaking "mass term" m r a 3 7 5 = m T /3r. The staggered chemical potential was 
already considered in jlO|, jll| to show the realization of irrational fermion charge induced 
by scalar fields and still maintaining the time-reversal symmetry. The parity-breaking mass 
matrix m T( 5r in this version of the graphene honeycomb structure is the continuous limit of 



the Haldane mass (energy) [12| linearized around the Dirac points jl3|- We consider it as 
constant as in jl^, but the same term was taken to be position-dependent (a domain wall) 



in 15j to investigate how the electronic properties of graphene are modified under such a 
domain wall. 

When such terms are added to the Hamiltonian, the corresponding Dirac operator an- 
ticommutes neither with a 3 nor with a 3 7s and the Hamiltonian does not admit any norm- 
preserving (conjugation) symmetry, that is, normalizable negative-energy eigenstates are no 
longer obtained from the normalizable positive-energy eigenstates and the lowest energy 
states are not necessarily self-conjugate under the operations of either a 3 or a 3 7s. Neverthe- 
less, it may be found a non-unitary operator which conjugates negative to positive-energy 
eigenstates and vice- versa, but, due to the lack of unitarity, the density of positive-energy 
states is different from the density of negative-energy states belonging to the continuous en- 
ergy spectrum and, as a consequence, an irrational fermion number may be induced. Since 
we do not know exactly the fermion single-particle energy states in a general magnetic field 
then we are going to consider the case of homogeneous magnetic fields whose single-particle 
energy levels are the familiar Landau levels. Once the energy spectrum is discrete, there is no 
way to realize irrational fermion number. Here the fermion number density is proportional 
to the surface density of states of the LLL whose associated eingenstates are eigenfunctions 
of either a 3 or a 3 7s, depending on the magnitude of the applied magnetic fields. 

Now, the Dirac equation is 

-iV - q^X - e~A) + + m T /3r]^(r) = E^(r), (22) 
10 



a.l 



where the vector potentials are written as in (j3J) and (1121) with A{r) and A§{r) given by 



A(r) 



Br 2 



and An(r) 



B,r 2 



with B, B 5 > 0. 



4 " w 4 

The set of equations for the spinor components can be written as 



(23) 



20* 

2a 



(E + mJ)il) a _, i 



2<a 



tp°L = (E — mJ)%l) b _ (24) 



where u± = eB±qB$ and m± = fi±m T . This problem is exactly solvable, but the spinorial 
structure of the eigenstates of the Hamiltonian operator depends on whether /i > m T or 
fi < m T as well as on whether eB > qB 5 or eB < qB 5 . Hereafter, we take fi, m T > 
and eB, qB 5 > 0, although the cases when eB, qB 5 < and/or fi, m T < can also be 
analyzed straightforwardly Even in the cases we are concerned here, there are four different 
situations, which can be better appreciated if examined separately. 

Case 1 : qB§ > eB . The energy spectrum is \E± \ = ^/2n \u±\ + m±, with n e N and the 
eigenstates are 



/ 



E + >0 



-i\E+\t 



V2 



\E + \ U " n > l > UJ + [ 





1 _ m + F? 

1 \E + \ rLn - 1 







E->0 



-i\E-\t 



V2 








(25) 



and 
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where 



E+<0 



E_<0 



J\E+\t 

~7T 



-,i\E-\t 



V2 



l-^R n ,i,„ + (r,9) \ 








Tra- 
ill 



■-n-n-i 



(26) 



Rn,l 




>±l 



/ i! Ve-^ r ^(|^ ± |r 2 /2) com I e N, (27) 

are the Gauss-Laguerre modes in cylindrical coordinates and expressed in terms of the 
associated Laguerre polynomials L^(|cu±| r 2 /2). Here R* stands for the conjugate to R whose 
normalization is J \Rn^ u ± (r, 9)\ rdrdO = 1, and we are assuming that R-i t i iU) ± (r, 9) = 0. 
For n = and for the case /i > m r (m± > 0) the lowest energy states are 



/ Ro,i,u + (r,9)\ 



ii+=m+ 



V 









/ 



E—=m- 



-im—t 



J 









\ 



(28) 



We notice that when m r — > 0, the energy states above are degenerate such that one has 
only one eigenspinor with the same structure found in 0171) . which is normalizable even when 
\i — > 0, that is 



( ip+(r,t) \ 





( Ro,i. 



V2 








V R 0,l,u 



(29) 



/ 



Had we started with m T = and \i < the LLL would have energy E Q = \x < and the 
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corresponding eigenspinor would be 








V2 



\ip b -(r,t)J 








(30) 



(r, (9) y 



Then the induced fermion density numbers, for [i — > 0, are given by: 



p=-(0|[^] |0) = -sgn(/i) 



47T 



and 



P5 = ^ (01^,75^] |0> 



-sgn(/x)— . 



(31) 



(32) 



The above results are compatible with those found in the previous section, f fl9|) and (}2T| , 
which reveals that the fermion numbers are proportional to the flux of the magnetic field. 
The dependence on the sgn(/i) is explained by the fact that when fi is taken initially as 
positive, the LLL is occupied by electrons on the conduction band (E = \x > 0), whereas 
when /i < the LLL is occupied by electrons on the valence band. The above result is only 
valid for \i — y 0, otherwise there would be a contribution (parity invariant contribution) 
coming from the higher Landau levels. 

Although there is no norm-preserving operator which conjugates negative to positive- 
energy eigenstates, when fi ^ 0, the energy spectrum is discrete and the eigenstates can be 
normalized, such that the density of states from the positive and negative energy spectrum 
(except for the isolated mode E = fi) is equal to each other. Then their contribution 
vanishes when one calculates the fermion numbers themselves and it is possible to show 
that the fermion numbers are proportional to the magnetic flux. In fact, the results in 
( 13T|) and ( |32|) express that the fermion density numbers are equal to the surface density of 
states of each Landau level (or the surface density of zero modes) and, as a consequence, 
the degeneracy is proportional to the magnetic flux. 

We emphasize that in this case, as was already revealed in equations f[T^]) and ( 12~T|) . the 
fermion charge is proportional to the flux of the (axial) magnetic field, whilst the fermion 
chiral charge is proportional to the flux of the magnetic field. This will be explored in the 
next section to analyse the effective action for the gauge fields in the relativistic context, 
where we will show that the derivative expansion approximation to the effective action leads 
to a crossed Chern-Simons term. 
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The next case we treat is well known and has been recently analyzed in by resorting 
to another representation for the reducible gamma matrices in 2+1 (space-time) dimensions. 
Although they do not consider the vector-axial gauge potential, their results are more general 
than ours in the sense that they take inhomogeneous magnetic fields and analyze also the 
behavior of the chiral condensates. 

Case 2: qB 5 < eB . The energy spectrum is identical to that in the previous case and 
the eingestates are expressed as 



E+>0 



E->0 



-i\E+\t 



r \E.. \ xh n,l,t»i+V i 



ijl 



\ 
( 



-i\E-\t 



Rn-l,l+l,u} + ( r , 









\ 



\E_\ n J>' 



y -Ul + |^|-Rn_i,j+i jW _ (r, 6) J 



(33) 



E + <0 
+ 



f\E+\t 



-Z A /l + ^Rn-l,l+l,aj+{r,< 



\ 
( 



£_<o e 



i\E-\t 



V2 








\ 



1 + ^Rn,i,u-{r,0) 



For fi = and m T > the eigenstates of the LLL is 

( ^ b + (f,t) \ 




E=rn T 







-im T t 



V2 



( R ,i, u+ (r,9)\ 








(34) 



(35) 
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and exhibits the same spinor structure as the one in (1161) . whilst for p = and m T < we 
have 



E=— \m T \ 





( ^ b + {f,t) \ 









e i\m T \t 










\ o ) 









V 



(36) 



/ 



Ro,i,u>- (r, 6) 


In obtaining the induced fermion charges it is convenient to separate the contribution 
of the LLL from the higher LL, which also contribute in this CcLSG clS well, but the last 
contributions vanishes when m T — > and one finds: 



p = -(o| |o> 



eB 

-sgn(m r ) — , 



P5 = (°l [^ f ,75^] |0) = -sgn(m T ) 



Air 



(37) 
(38) 



which are consistent with (II S p and (12U|) . 



IV. INDUCED FERMION CHARGES FROM FIELD THEORETICAL CALCULA- 
TIONS 

As one knows QED in 2+1 space-time dimensions can be enriched by adding a topological 



mass to the gauge field 16] and that such a C-S term can be generated by means of per- 
turbation calculation jl7|. The C-S term comes from the first-order in external momentum 
contribution of the vacuum polarization diagram. 

Here, we show that C-S terms can also be induced in the extended QED, whose La- 
grangian density is 

C = * [r(d, + eA, + gAf) - p^ - m T r]^. (39) 

We also compute the fermion induced currents from the C-S effective action. 

The vacuum polarization diagrams we have to consider are those shown in FIG. [T] which 
correspond to the vacuum polarization operator 



Dr(*) = -i I |^tr 



— 7 M (e + g7 5 )-7 3 Y(e + g 7 5 ) 



(40) 

where Sp(p) = {j/> — p^3 — m r r)~ 1 is the fermion propagator given in terms of momenta and 
under a staggered chemical potential and parity-breaking mass. 
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FIG. 1: 



Here, as in the previous section, we are going to consider the contributions of the staggered 
chemical potential and of the parity-breaking mass separately. This is not only because we 
want to see the dependence of the fermion currents on the signs of [i and m T as in the 
previous section, but also because the fermion propagator can be easily handled when the 
approximations 

S F (p) « S F>ll (p) = , 1 3 » m T ) (41) 

f-vn 

and 

S F (p) ~ S FtmT (p) = - — (/x « m T ) (42) 

p — m T r 

are employed. 

In both cases we look for the term which is first-order in the external momentum, namely 
n^" ~ e ixva k a . By taking the first approximation for the fermion propagator, Eq. (|4ip . we 
note that the third diagram in FIG. [T] is the only one which gives a contribution to the C-S 
terms. This can be appreciated through the following expression for IP" 

and by noting that the Levi-Civita tensor comes from tr(7 At 7 1/ 7 a r) = — Aie^ ua . By isolating 
the relevant terms to the C-S contribution one finds 

= -^e^ a k a \im [ dz [ d 3 p— ] 2 = 

27T 3 fc->o J J [p 2 - fi 2 + k 2 z(l - z)] 2 

= i^le^k hm r - = i^-eT'k ± (44) 

From the effective action up to second-order on the gauge fields corresponding to the 
third diagram 

S 3f = ^J d3x {M x ) nia, ( k * = ^ a )A^(x) + Af{x)W v (k a = id a )A v {x)} , (45) 



(p + k) 2 — fi 2 



7"(e + g7 5 )- 



p z — ji 



6 

2 r 



e + q~/ 5 ) 



(43) 
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one obtains a crossed C-S term, namely 

Si 2 f f = e^ a I d 3 x {A,(x)d a Af\x) + Af{x)d a A v {x)} . (46) 

I I J 

The induced fermion current density is obtained from this contribution to the effective 
action through 

if) = -rf— = -f-^r" 9 ^ 

oA^ Air \fj,\ 

- u 7) 

- 4tt i/xr (5) ' 10 

where = {—l/2)e^ ua F ya ^ is the dual of the (axial) field strength, whereas the (axial) 
fermion current density can be obtained as 

= ( 48 ) 
Air \fi\ 

By taking the time-component of (l4"Tj) and ( 148]) one recovers the results for the fermion 
charge densities, equations (l3Tj) and f )32|) . respectively. 

Now, we consider the effective C-S action coming from the contribution of the first two 
diagrams. We notice that such contributions are obtained by taking the second approximate 
fermion propagator given in ( I42p . Here we calculate both contributions together, that is 



n^(fc) = -i I T^tr 



(2 



7T 



— ^ -7 M (e + g 7 5 )^ - 7 " e + g 7 5 ) 



(49) 



(p + k) 2 — m 2 p 2 — m 2 

This time there will be no crossed C-S terms because tr( 7 M 7 v 7 a T) = —Aie^ ua is obtained 
from the contributions proportional to m T and to e 2 and q 2 , thus one has 

Kst = «'i^T^^ a *« ^ fi r 2 n, = (50) 

for the first and second diagram, respectively, and the corresponding C-S effective action is 
~Sfff = j^| ^ Q J d 3 x {e 2 A,{x)d a A u {x) + q 2 Af{x)d a Af{x)} . (51) 
From (lo"Tj) one obtains the induced fermion currents 

n*\ - 6 J^1 - -^^ F ^ an d m - 5 -^L- - -^J^ F » 
{J) ~ SA„ ~ 47r|m T f (J5) " 5 4 5) " 4vr|m T r( 5 )' 

whose time-components recovers the results (1571) and (J5S|) respectively. 
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V. CONCLUSIONS 



We have obtained induced fermion charges in a continuum chiral theory for massless 
planar electrons in a honeycomb structure by means of only vector gauge potentials. For 
some configurations of the vector gauge potentials, particularly quantized magnetic vortices, 
we show that the induced charges can have fractional values. Since isolated and localized 
zero-energy states together with symmetry of sublattices are crucial for the appearance 
of fractional charges, we have analyzed the spinorial structure of the zero-modes as well as 
the spectral-symmetry conjugation of the first-quantized Hamiltonian. Moreover, in order to 
understand and highlight the origin of zero-energy states, that is, if they come from electrons 
on either valence or conduction bands we have studied in detail the energy eigenstates of 
fermions in homogeneous magnetic fields by taking also into account a staggered chemical 
potential (/x) and a parity-breaking mass term (m T ), which open a mass gap in the energy 
spectrum (Landau levels). By taking /i,m T — > 0, we calculate the induced fermion charges 
and the results are in consonance with the results found when the electrons are under the 
influence of finite magnetic fluxes. 

In pursuit of the realization of quantum anomalies in condensed matter systems, we also 
also discuss, through field-theoretical calculation, the relation of induced fermion currents 
with the parity-anomaly in an extended Quantum Electrodynamics which involves a vector 
and a vector-axial gauge fields. Unfortunately, this relation is not established beyond doubt, 
since the parity-anomaly is only realized if parity symmetry-breaking terms, such as the 
staggered chemical potential and the Haldane mass, are present in the physical system, and 
that is not the case in graphene. 
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